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Abstract 
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corrections are taken into account, new effects in neutrino oscillations, as — > i/^, appear. This 
oscillation is constrained by the solar neutrino data. In our analysis, we have found an allowed 
region for our two free parameters e and m\. The radiative correction, e, can vary approximately 
from 5 X 10~^ to 10~^ and the calculated fourth mass eigenstate, m4, 0.01 eV to 0.2 eV at 2a level. 
These results are very similar to the ones presented in the literature. 
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I. INTRODUCTION 



Over time, we know more about the parameters of neutrino oscillations and such knowl- 
edge is crucial for us to restructure the Standard Model. For a recent statistical analysis 
of all experimental neutrino data available see [1, 2]. Moreover, since the LEP data, we 
know that there are only three active neutrinos [3]. However, from several experimental 
results and theoretical speculation, there would be a new kind of neutrino that must be not 
sensitive to the interactions of the Standard Model, i.e., it would not couple nor to the Z 
and neither to the W . They only couple, in the context of the known physics, to gravity, 
and for these reason they are traditionaly called "sterile neutrinos" . For a recent review see 
Refs. [4]. On the other hand, these sort of neutrinos might interact to non-SM particles. 
Once sterile neutrinos are added they can be of several types depending on the mass scale 
related with them. They may be or not related to some anomalies [5] in neutrino data [6-9] 
or with the results of WMAP-7 [10], which indicates the existence of four relativistic species 

If sterile (right-) handed neutrinos are introduced, an interesting possibility is that a Majo- 
rana active neutrino is mass degenerated with a sterile Majorana neutrino, since in this case 
the two Majonana neutrinos are equivalent to a Dirac massive one. When quantum correc- 
tions are considered the mass degeneracy is broken and we have a gnasi-Dirac neutrino [11]. 
We distinguish this sort of neutrino from the ones called i^settdo-Dirac neutrinos [12]. In 
the latter case, the mass degenerated neutrinos are both active. It is well known that solar 
neutrino data impose strong constraints on the allowed mass for this sort of quasi-Dirac 
neutrinos [13]. Here we shall apply these solar neutrino data for the particular model of 
Refs. [14, 15], which is an interesting implementation of the model of Refs. [16, 17]. In 
these models the PMNS mixing matrix in the lepton sector is given, at the tree level, by 
the tribimaximal mixing matrix. Also, these models have enough interactions to generate 
corrections to both masses and mixing angles that we can expect a new channel for the 
disappearance of electron active (anti-)neutrinos. Moreover, these corrections may be, or 
not, responsible for the deviations from the tribimaximal mixing matrix [18] in the leptonic 
sector as some experimental results indicate [19-21]. 

To implement this program, we will apply the available solar neutrino data to the fol- 
lowing channel: electronic neutrino, Vf.^ oscillating into sterile neutrinos, say Vg. For the 
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statistical analysis, our model has two parameters: one mass eigenstate (mi) and the radia- 
tive correction (e). This is an important difference if we compare our analysis with the one 
made by de Gouvea et al. in section (III. A 2+1 case) [13]. Despite the differences in the 
model building, we have obtained an allowed region for e that is similar to the corresponding 
values found in [13] . Also important to point out that we want to constraint our quasi-Dirac 
model by the solar neutrino data and not to solve the experimental anomalies presented 
in [6-9] , however our mass scale is consistent with WMAP-7 [10] results. 

The content of this article is the following: first, in Sec. II we review the basic features 
of our quasi-Dirac model, showing its basic structure and possible interactions with the 
correspondent radiative corrections to the neutrino mass matrix. In Sec. Ill we briefly 
describe all the solar neutrino experiments and their respective results. Also we include 
the method of statistical analysis used. In Sec. IV we show in what condition the model 
is consistent with solar data, showing the allowed regions in the parameter space (mi,e) 
for our quasi-Dirac model and discussing the results. Finally, some concluding remarks are 
present in the last Sec. V. 

II. THE QUASI-DIRAC SCHEME 

Recently it was shown that it is possible that all neutrino flavors are part Dirac and part 
Majorana [14, 16]. The latter occurs because two of the four Majorana neutrinos are mass 
degenerated and have opposite parity, so they are equivalent to one Dirac neutrino. The 
other two have distinct Majorana masses. It is also well defined in the literature, as we 
said before, that when these two neutrinos form a Dirac state and they are active, we called 
them pseudo-Dirac neutrinos. When there are one active and one sterile, they are called 
quasi-Dirac neutrinos. In our particular model, we point out that there are initially three 
right-handed neutrinos. Two of them are integrated and we obtain a model "3 -|- 1" - three 
active neutrinos and one sterile [22] . 

In this section we are briefly going to describe the construction of our model (Sec. II A). 
In subsection II B, we show the main interactions that are going to be used to obtain the 
radiative corrections for the neutrino masses. These radiative corrections are very important 
to our analysis: we study their effects on the break of the degeneracy between the two 
Majorana neutrinos that form a Dirac neutrino at tree level. 
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A. Model Building and Basic Elements 



The model we are going to present here is based on a gauged B — L symmetry with a 
quasi-Dirac neutrino in which the right-handed neutrinos carries exotic local B — L charges. 

When the ^'3 symmetry is added to the model, the left-handed leptons belong to the 
reducible triplet representation (3 = {L^, L^, Lr)) since all of them have the same B — L 
charge. The triplet can be decomposed into irreducible representations as 3 = 1 + 2, then 
we can write the singlet and doublet as follows: 

L2 = -^(Le + L^ + L,) - 1, 

(Li, L3) = (^(2^e -L,- L,), ^(L, - L,)^ ~ 2. (1) 

On the other hand the right-handed neutrinos have different B — L charge, so they can 
transform under S3 only as a singlet 1 = with B—L — —4, and a doublet, 2 = {rieR, tItr), 
with B — L — 5. The scalar sector has two scalar doublets of SU{2) with weak hypercharge 
Y — —1 that are denoted by $1,2 = (<^S^2 <^i^2)^- They are singlets of and we will denote 
((/9°((/92)) = ^1(^2) / \/2. If n^R is considered light but UeR and Utr heavy, we will integrate out 
the latter ones. After that the effective lepton Yukawa interactions are given by a dimension 
five effective Lagrangian plus a Dirac mass term as follows 



-Cf = h,L,^,n,R + ^ [(^li)r%][L,l%] + ^ [{LI)r%][L3l%] + H.c, (2) 

where the mixing angles in the {ueR-i Utr) sector have been absorbed in the dimensionless 
couplings, /i2 and h^. 

Prom the Yukawa interactions in (2), we obtain the mass matrix in an appropriate basis. 
At the tree level the mass matrix is diagonalized by the following 4x4 matrix: 

Uo=i ^™ \ , (3) 

\ Oix3 1 / 

where Utbm is the tribimaximal matrix and denote the matrix row or column with entries 
equal to zero. 

Since the model has more interactions than those in the Standard Model, the neutrino 
mass matrix, when radiative corrections are taken into account, is not necessarily to be 
diagonahzed by Uq, written in Eq. (3), but for Uq — >■ U, where U is now another 4x4 
matrix. We will not write this explicitly, because it is not the aim of this manuscript. 



B. Quantum corrections 




FIG. 1: A 1-loop contribution for the Major ana mass matrix in the flavor basis induced by charged 
scalars. The • in the left vertex denotes interactions in Eq. (5). 



When radiative corrections are taken into account - see Fig. 1 - the neutrino mass matrix 
can be written as 

M" = M^" + AM^ (4) 

where M^'^ is the mass matrix at tree level [15] and AM*^ arises from 1-loop corrections. In 
order to calculate the mass corrections, AM^, we have to consider all the Yukawa interactions 
in the lepton and scalar sectors. From Eq. (2) the scalar-charged lepton interactions are: 

1 mi 



-^CC = ^^i^L + f^L + TL){iyD)R(l)t - ^ 



i[2(e-)^ + (/i^)^ + (r^)^]z/iL 

1 



^{v; + Re0°* + ^Im0°*)0+ + H.c, (5) 



and we have used 1^^3)1 = i^{i,3)l, ^2L = ^21, a/2(z/^)/j = n^R. 

The Yukawa interactions from which the charged leptons get mass are mainly the diagonal 
ones, 

^yuka^a ^ G J^i^l 1){^iIr + H.C, (6) 
I 

where = and / = e,fi,T. The charged lepton mass matrix is almost diago- 

nal [14], hence m/ ^ Gvi, and the neutrino interactions with charged leptons are given by 
{mi/vi)uiLlR4>^ ■ We stress the fact that it is the neutrino flavor basis that is important here. 
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On the other hand, the scalar potential includes the following interactions: 

$2, • • • ) « Xii^l^M^i + Aa$}*2*^*z + K^^{e^sM<Px + Ky^e^.^c^y + H.c, (7) 



where I = e, ji, r. In Eq. (7), ^^^^ denotes a scalar doublet with Y = +1 and without B — L 
charge, and (p^i 4>y are scalars carrying also B — L charges [14]. 

With the interactions in Eqs. (5), (6) and (7), we obtain diagrams hke the one in Fig. 1. 
As we mentioned before, these sort of diagrams provide corrections to the Majorana masses 



for the active neutrinos, i.e., iyaiYi^hL- Also corrections to the Dirac mass terms VaLn^R arise 
from diagrams similar to the one shown in Fig. 1. The 1-loop corrections to the neutrino 
mass matrix [M^) in = N'^j^ + {N[^Y basis where N'^^ = {i/f.^n^r^'j^)^ are written in the 
following: 



I 1(1 -^e) 



-^(l-2ee) 



-1(1 -2e, 



(1 + <) \ 



-^(l-2e,) 
-1(1 -e,) 



'1 _|_ ms 
^6 ~'~ 2m 



2m 



, fl + eM 



From Eq. (6), G = rrie/ve = m^/vn = rrir/vr, and we can express 

^-i-^^Ai,m,^ ABi[mni - ln(m,VmJ,)] - mj, + Q ln(G/mJ,)], 



(8) 



^^WAS.fmni - Hml/ml,)] - + ain(G/m^,)]. 



(9) 



When all e's and e"s in Eq. (8) are equal to zero, the mass matrix is the same as the one 
represented at tree level. 

In Eq. (9), A, Bi and Ci are given by: 



^ ^ V^hky{(l)^){(l)y) ^ 



V2 



rrr 



2(mf — m? + )('mf — rri^ + )imf — m'^ + )(mf — m?o ) 



Ci — m + +m + +m + +m + — brrii . 



(10) 



The general form of the mass matrix in Eq. (8) is very complicate to treat, so we will do 
some approximations in order to simplify our analysis. We consider two cases: 



1. CASE A: Aiem! — Ai„m?, = X^-^ml = and Xy.^m^ — X9.,,ml — A^^m: 



2e' 



^2^"V 



\2t' 



m 
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In this case we have eg = = = e and e'g — e'^ — e'^ = e' where 

_ \/3 vi ky{(j)^){(l)y)m? 

e'^%e. (11) 

n 

M is a typical mass in the charged scalar sector and m^o is the mass square of the 
Higgs of the SM. 

2. CASE B: Xie^e = Ai^m^ = XirTUr = m. 

In this case we have a certain hierarchy in the radiative corrections (eg < < e,-), 
represented by: 

_ me _ m^ 

4 = ?4 < = ?4 (12) 



where 



1 vi 

mrrirA, 



8V^7r2 V2 m\ 

VSM 

1 V2 M 



2 



e' ~ 5 — mrrirA. (13) 

We are going to use both of these approximations in the following analysis. Therefore, 
we have two main free parameters which were written in Eq. (8): the mass mi (in eV units) 
and the radiative corrections (adimensionally), e for the CASE A; and €-,- — € and e'^ — e' for 
the CASE B. For this case, the other e's and e"s are calculated by Eq. (12). We notice that 
e-,- fa e'^, then we are going to express our results, for the CASE B, using e^- — e. 

For more details about this quasi-Dirac model building, we suggest the reader to check 
references [14, 15]. 



III. SOLAR NEUTRINOS CONSTRAINTS 



The detection of neutrinos travelling from the sun has given us a tremendous evidence 
of neutrino oscillation. We might say that it was the first time that physicists were doing 
astronomy with neutrinos and several aspects of the solar behavior have being observed and 
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understood since then. Prom Homestake to SNO, nowadays we have an amount of consider- 
able and significant data that also gives us the opportunity to use this fact to constrain and 
test the validity of models. This is exactly what we are going to do: constraining the param- 
eters of the quasi-Dirac model presented in Sec. II and checking its validity in confrontation 
with the solar neutrino data. This data is taken from the following experiments: Homes- 
take [23], Gallex/GNO [24], Sage [25], Kamiokande [26], Super-Kamiokande [27], SNO [28] 
and Borexino [29] . In Sec. Ill A we are going to present a small review about these experi- 
ments and their main numerical results. In Sec. Ill B we present how to treat the oscillation 
physics of solar neutrinos and the main points of our statistical analysis. For recent reviews 
on solar neutrinos see [30-32]. 

A. Experimental data 

For the statistical analysis. Sec. IIIB, we are going to consider the entire set of the solar 
neutrino data presented in Table I. This table presents each solar neutrino experiment and 
the measured fiux {(pexp)- Depending on the experiment, the fiux is measured by charged, 
neutral current reaction and elastic scattering. 

Elastic scattering experiments, z/^ + e^ — > z/^ + e" (a = e, /i,r), include Kamiokande [26], 
Super-Kamiokande [27], Borexino [29, 34]. Kamiokande and Supcr-Kamiokande detected 
^B neutrinos with threshold of 7.5 MeV and 5 MeV, respectively. Borexino, on the other 
hand, detects neutrinos from the ^Be fine with an energy of 0.86 MeV. Recently, Borexino 
has also measured for the first time the fiux of low energy pep neutrinos: 0pep = (1.6 ±0.3) x 
10^ cm~^s~^ [35]. However we do not use this values in our analysis. 

SNO experiment [28] detects electronic neutrinos in a charged current reaction, z/g + (i — )■ 
p + p + e~ (threshold of 5 MeV). Also, there is the detection of other neutrino flavors by 
neutral current reaction, Ua + d^n + p + i/a (threshold of 2.225 MeV), and elastic cross 
section. SNO had three stages and obtained different fluxes [36-38] shown in Table I. 

We used also the Homestake experiment [23], Ug+^'^Cl — > ^^Ar -|-e~ (threshold of 
0.814 MeV), and the ^^Ga experiments: GALLEX/GNO [24] and SAGE [25] (threshold 
of 0.233 McV). In Table I, we referred to all gallium experimental results [39]. We notice 
that they are sensitive to almost the entire neutrino solar spectrum. 



8 



Experiment 


Experimental Data 


Homestake [23] 


2.56 ±0.16 ±0.16 SNU [33] 


Gallex/GNO and Sage [39] 


68.1 ± 3.75 SNU 


Kamiokande [26] 


(I^Kam = (2.80 ± 0.19 ± 0.33) X 10^ cm-2s-i 


SK [27] 

L J 


(l)SK = (2.35 ± 0.02 ± 0.08) x 10^ cm-^s-i 


SNO - D2O [36] 


4>cc = {l.76t'omi^tat.)tr,l{syst.)) x 10^ cm-^s"! 
(f)ES = {2.39l:oitistat.)toil{syst.)) x 10^ cm-^g-i 
(bNC = (5m+°r,it(stat.)+^AUsyst.)) x 10^ cm-^s"! 


SNO - NaCl [37] 


(t>cc = (1.59tn?(^*«*-)lo:o8(«y^«-)) X 106 cm-^s-i 
= {2.2ltoil{stat.) ± 0.10(syst.)) x 10^ cm-^g-i 
(bNC = (5.21 ± 0.27(stat.) ± 0.38(syst.)) x 10^ cm-^g-i 


SNO - ^He [38] 


<^CC = {l.67t''o^olistat.)tZlisyst.)) x 10^ cm-^-^ 
= (1.77tOit(stai.)iHo(«y^«-)) X 106 cm-2s-i 
= (5.54l0:33(stat.)l°J6(syst.)) X 106 cm-2s-i 


Borexino [29] 


= (4.84 ± 0.24) X 10^ cm-^g-i 



TABLE I: Resume of solar neutrino data. Also including the experimental uncertainties. 1 
SNU=10~^6 captures/atom/sec. 

B. Analysis 

Neutrinos are produced for several thermal nuclear reactions in the center of the sun [40] 
and we present them in Table II extracted from [41]. The energy of these neutrinos are 
of a few MeV. To be accurate, we must treat the center of the sun as a region where the 
chemical composition modifies itself with the radius. So each reaction produces a different 
flux of neutrinos and this changes, as we pointed out, with the position from the center 
of the sun. Neutrino sources are the so-called: pp, hep, pep, ^^N , ^^O, ^'^ F , and ^ Be. 
Details of the distribution of the neutrino production as a function of the radius for each 
of the solar neutrino sources can be found in [40] and we used this profile in our work to 
average the oscillation probabilities, since detectors only "see" these averages. 

After electronic neutrinos (fe) are produced by several reactions and in different points 
of the core of the sun, they will propagate inside the sun, which has a radius Rsun ~ 
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Source 


Reaction 


Average v Energy (McV) 


Maximum v Energy (MeV) 


pp 


p + p^d + e'^ + Ve 


0.27 


0.42 


pep 


p + e~+p^d + Ue 


1.44 


1.44 


hep 


^He ^-p^"^ He + e+ + Ve 


9.63 


18.78 


^Be 


e- Be -^'^ Li + 


0.86 


0.86 


8B 


Be* + e+ + Ue 


6.74 


15.00 


13n 


13^ ^13 C + e+ + Ve 


0.71 


1.19 


150 




0.99 


1.73 


17p 


17^ ^17 + e+ + Ve 


0.99 


1.74 



TABLE II: Sources of solar neutrinos: first column represents the name of the source which pro- 
duces the electronic neutrino inside the sun; the second column shows the resume reaction; the 
third and fourth columns represent, respectivelly, the average neutrino energy and the maximum 
neutrino energy. 

6.9 X 10^° cm. This propagation is described by the effective Hamiltonian of the system in 
the flavor state base: 



^e//(r) 



IE 



(14) 



We emphasyse that — M''{e,mi) is taken from Eq. (8), V{r) is the potential of the 
neutrino interaction with the solar environment, E is the i/g energy and r is the distance 
from the center of the sun. The potential, V{r), can be written as the sum of the charged 
current and neutral current interaction {V{r) = Vcc{r) + Vnc{r)), which are dependent of the 
electronic density (ne(r)) and neutron density (n„(r)) of the environment. Both of these 
quantities change with the distance from the solar core and can be written as: 



V{r) = Kc(r) + Kc(r) = V2GFn,{r) - ^\/2Gfn„(r). 



(15) 



The profile of ne(r) and n„(r) that was used in our analysis has been extracted from [40] 
and Gp is the Fermi coupling constant. 

The survival probability {Pee), for each energy and in each point of neutrino production. 
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is calculated from the amplitude Age, which can be written as: 



X 



( 1 ) U^ac X diag{exp (-i^'J, exp {-i<^>'^) , exp (-i^'g), exp (-^$4)) 

/ 1 n n \ 

x/r^ 







Pc -VPc 



X (iia5f(exp(-i$i),exp(-i$2),exp(-i$3),exp(-i$4)) x t4 



mat 








(16) 



So the survival probability is written as Pee = l^eeP- In Eq. (16), Umat = t^mat(e, mi, £■) 
is the matter mixing matrix which diagonalizes the effective Hamiltonian represented by 
Eq. (14). The crossing probability, which will be discussed later, is represented by Pc- The 
matrix U^ac = U^aci^^f^i) is the vaccum mixing matrix which diagonalizes Eq. (14) when 
V{r) is equal to zero (vacuum regime). If we take e = 0, no radiative corrections, for any 
mi value, U^ac is going to be the tribimaximal mixing matrix. As the elements of U^ac, the 
elements of Umat are modified by the choice of the parameters e and mi. When the electronic 
neutrinos travel to less dense regions of the sun, Umat U^ac- The phases represent the 
evolution of the mass eigenstates in matter. We express this as: $i = J^^"""" ijf{x)/{2E)dx, 
where fJ^f{x) {i = 1,2,3,4) is the mass eigenvalue of Eq. (14) and Tq is the neutrino point 
of production. The phase <l>'j has a similar meaning to $j, but for the vaccum propagation. 
We discuss about it later in this section. 

In Fig. 2, we show the evolution of the mass eigenstates in the sun for a neutrino with 
energy E — 5 MeV, mi = 0.001 eV, e — 1.0 x 10~^. The sohd black curve represents the 
mass eigenstate /^i; dotted blue, dashed green and dotdashed red ones represent 112, f^4 and 
/X3, respectively [42]. We notice that U2 and 2/4 are practically degenerate, which is the most 
important characteristic of quasi-Dirac models. It is also possible to notice that matter 
can break this degeneracy for very small radius as we can notice in Fig. 2. However, for 
very small e, we see that 1/2 and 1/4 are practicaly degenerate, generating a Am^a that can 
be sensible to oscillations: Aml2L/{2E) ~ 1, where L is the Sun-Earth distance, which is 
about 150 million kilometers. 
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fif /2E (eV) 




FIG. 2: Evolution of the mass eigenstates inside the Sun. For this plot we use e = 1.0 x 10~ , 
E = 5 MeV and mi = 0.001 eV. Black solid curve is for ni; blue dotted one for fi2; green dashed 
curve for /i4 and red dotdashed is for ^s- 

For instance, in the limit of e — > 0, we recover the original and standard 3x3 situation 
without the sterile neutrino presence, where the terms f/ei and Ue2 solve properly the solar 
neutrino problem: the deficit of z/g arriving the Earth. For e 7^ and small, it is important 
to notice that z/2 and will be a coherent mixture - (z/2 + iu^) / - and this is an eigenstate 
of the Hamiltonian in vacuum. 

We know that the mass eigenstates can feel MSW resonances during the propagation [43] . 
When neutrinos go through the MSW resonance, the conversion probability is maximal. We 
remember that conversion probabilities are obtained using the expression written in Eq. (16), 
but changing the position of the number "1" of the line vector (1 0). For example, Pe/^ is 
obtained using the line vector (0100). In principle, we can have resonances among all the 
mass eigenstates, however, z/3 is the heaviest and it will not suffer resonance - its propagation 
is adiabatic. Also, we can say that the scale Am|j, with i = 1,2,4, can be averaged out. 
In other words, this mass squared difference scale is not important for the solar oscillation 
phenomenon. The moment of the resonance is represented by the matrix that contains 
Pc in Eq. (16). This Pc is the crossing probability, which represents the probability of a 
mass eigenstate z/j be converted to another mass eigenstate Uj. In the standard neutrino 
oscillation case, if the propagation is adiabatic, we must have Pc = 0. In other words, there 
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is no conversion between two mass eigenstates. In the instantaneous mass basis (i/f*, for 
i — 1, 2, 3, 4), where is the neutrino state in matter, the evolution equation is expressed 
as: 

fli/^ r 1 _ _ _ _ r^nJ-rM 

(17) 



dx 



where n^{x) is the effective mass eigenstates calculated from the eigenvalues of Eq. (14). If 
the last term of Eq. (17) is significant compared with the first one, non-adiabatic transition 
can happen. The adiabaticity parameter, represented by the letter 7, is evaluated at the 
resonance point, for simplicity, and is defined as: 



(M|(a:)-/.f(x)) 
2E 



[Ul{x 



dx m 



(18) 



When << 1 {jij >> 1), the propagation is non-adiabatic (adiabatic). Considering i/i and 
1/2 (or 1/4, since they are practicaly degenerate), for any values of e and mi, and evaluating 
Eq. (18), we conclude that MSW resonance and all the propagation is adiabatic. Actually, 
for e — > 0, mixing angles and mass squared differences extracted from our model are very 
near to the experimental ones [3], so we know from experiments that the propagation is 
adiabatic. As an example, we can see in Fig. 3, evaluated for = 5 MeV, e = 0.8 and 
mi = 0.006 eV, that 712 is very large and much greater than 1. Then we can say that Fig. 3 
has showed, even for large e, that the value of 712 is kept large and then the propagation 
remains adiabatic. Even for larger values of mi we obtain the same pattern and magnitude 
of 7i2. We also can say that Fig. 3 represents the behaviour and magnitude for 714. However, 
since we have (quasi-) degenerate state between 1/2 and 1/4, 1712 ~ m4, we cannot say that 
this transition is always adiabatic. So, 1/2 ^ 1^4 is very dependent of the values of e and 
mi. We need to compute Eq. (18) and calculate the crossing probability, Pc, for this kind 
of transition. 

The crossing probability can be written as [44] 

p-l\U(vac)2i? _ p-1 

Pc = —, (19) 

where C/(^ac) 24 is the 24 element of the mixing matrix in vacuum and 7 = 724. If 7 is large, 
we have an adiabatic propagation of 1/2 and 1/4, and they will get out independently of the 
sun, as distinct mass eigenstates. Then, in this situation, Pc — 0. On the other hand, with a 
very small 7, which generally happens for very small e, we have a non-adiabatic propagation 
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FIG. 3: Parameter of adiabaticity for the ui — > 1^2 transitions, using e = 0.8, = 5 MeV and 
mi = 0.006 eV. 

of 1/2 and 1/4^, and they will get out of the sun, as mentioned before, as an coherent mixture. 
Then in this situation, Pc = 0.5. 

After the propagation inside the sun, neutrinos will travel in vaccum, with a phase 
i = 1, 2, 3, 4 - see Eq. (16) - where ^[ = J / {2E)dx and rrii is the neutrino mass eigenstate 
in vaccum. The integration is taken along all the path from sun to Earth. Only Am|4 is 
considered, for certain values of e and mi, and cannot be averaged out. Other mass squared 
differences are averaged out since Am'^L/{2E) >> 1. It is important to stress out that 
we are ignoring Earth matter effects. An extended analysis involving four neutrino families 
with a great variety of Am^^ and mixing angles was done in [45]. 

The general expression of the expected event rate in the presence of oscillations in exper- 
iment j in the four neutrino framework is given by Rf^: 

= 5^ '/•fc / dE\k{E) X [ae,i(E)(Pee) + (T,,,(E)(1 - (Pee) " (Pea))], (20) 
fc=l,8 

where E is the neutrino energy, 0^ and are, respectively, the total neutrino flux and 
the neutrino energy spectrum normalized to one from the solar nuclear reaction k with 
normalization given by the model BS05(OP) in [40] - see Table II. In Eq. (20), ae,i {cFx,i) 
is the Vf. {vx-iX — II, t) interaction cross section in the Standard Model with the target 
corresponding to experiment j, {Pee) is the average survival probability in the production 
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point, {Pea) and {Pes) are, respectively, the average conversion probability in the production 
point of Ve {O' — 1-1', t) and — >■ i/g. 

The test is calculated by: 



where Ri^p is the experimental rate for j-experiment - see Table I. Generally, the rate is 
defined as R = <p/<pssNh where <pssM is the total flux of the solar standard model extracted 
from the model BS05(OP) [40]: Rth corresponds to a (pth that represents the oscillated 
flux, which is related to the parameters of our model and evaluated using Eq. (20); Rexp is 
based on the flux 0exp, which is the experimental value extracted from Table I. The is 
calculated for each set of e and mi. Note in Eq. (21) that crj is the error, which takes into 
account the experimental error of a particular experiment and errors associated with the 
flux expectations in BS05(OP). 

IV. RESULTS 

Prom Eq. (16), we plot the probabilities for two sets of parameters (mi, e) for the CASE A. 
In Fig. 4, we used as input e = 5.0 x 10~^ and mi = 0.003 eV. The masses m2 and ma, which 
appear in Eq. (8), are calculated. To maintain hierarchy, which we choose the normal one 
for simplicity, the masses m2 and ma will be written as following: m2 — \/ Am^^„ + and 
ma = y^Amlt^ + mj. The values of Am^„„ = 7.58 x 10"^ eV^ and Am^^^ = 2.35 x lO'^ eV^ 
are the best-fit values at la taken from [3]. 

After introduce these values of masses in Eq. (8), we evaluate the neutrino evolution 
Hamiltonian in matter, Eq. (14), calculating the new mass eigenvalues and diagonalizing it 
to obtain the new mixing matrix. It is important to notice that for each set (e,mi) we are 
going to have new elements of the mixing matrix and new mass eigenstates etheir in vacuum 
and in matter. Then we evaluate the probabilities using Eq. (16) and average them in the 
region of production. In Pig. 4, the survival probability, P^e, is represented by the solid 
black curve. Conversion probabilities are represented in the following way: Pea {a — A*, t) is 
the dotted red curve and Pes is the dashed blue curve. 

In Pig. 5 we have plotted the probabilities for the CASE A for e = 1.0 x 10~^ and 
mi = 0.003 eV. Notation and representation of the curves are the same as in Pig. 4. 




(21) 
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FIG. 4: Survival probabilities (Pee - solid black curve) and conversion probabilities {Pea - dotted 
red curve - and Pes - dashed blue curve) for the CASE A. All the curves correspond to e = 5.0 x 10^''' 
and nil = 0.003 eV. We do not show these probability curves for the CASE B, since their behaviour 
is very similar and there are very small differences. 
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FIG. 5: Survival probabilities (Pee - solid black curve) and conversion probabilities {Pea - dotted 
red curve - and Pes - dashed blue curve) for the CASE A. All the curves correspond to e = 1.0 x 10^^ 
and mi = 0.003 eV. Similar results are obtained in CASE B. 

We notice, for example, for small e, such as e = 1.0 x 10~^ in Fig. 5, that the quasi-Dirac 
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situation mimics the standard one, since we do not see a significant conversion to the sterile 
neutrino fiavor. So, when e — > 0, or simply to very small values, we approach the traditional 
solar neutrino solution. However, for higher e, such as we saw in Fig. 4 (e = 1 x 10~^), 
conversion to sterile neutrinos can be significant for the entire neutrino spectrum. If we get 
an even higher e, we will see an even larger oscillation pattern of Pes- This also happens 
for the other channels of oscillations (fg — >■ i'ij.,t)- That is because with a larger radiative 
correction, e, we get Am'fjL/{2E) >> 1. The phenomenological effect of mi is very similar. 
For large values of mi, if we fix e 7^ 0, we are going to have more oscillation if we compare 
to the situation with a smaller mi. 

One of the main sources of neutrinos, considering SNO and SK as experiments, is the '^B. 
For energies above a few MeV, SNO and SK reveal that Pee ~ 0.3 and also Pea ~ 0.7. This is 
a very strong constraint. For energies below 1 MeV or so, the constraints come mainly from 
Borexino, Homestake and the gallium experiments. Borexino imposes Pge ~ 0.51. Then if 
we have larger values of Pes, Pee must be higher to compensate the disapeareance of active 
neutrinos {v^ or v^.) that would arrive in Earth detectors. For even lower neutrino energies, 
mainly of the pp chain, gallium experiments impose Pee ~ 0.5. So any modification on Pes 
for the set of parameters e,mi has to be compensated by Pee, especially in the high energy 
part of the spectrum. 

We do not show the plots for the CASE B, because the behaviour and the pattern of the 
curves are very similar. 

Next, we proceed a fit to the data. This will be used to constrain our model for 
both CASES A and B. We can define Ax^ = - = 0), where = 0) is 

valid for any value of mi, since e = represents the standard situation and solar neutrino 
experiments are sensitive only to the mass squared difference and not to the absolute value 
of neutrino masses. For e 7^ 0, we choose mi to vary from 0.001 eV to 1 eV. We remember 
that Katrin will impose a superior limit on neutrino mass of about 0.2 eV [46]. 

In Fig. 6, we present the allowed region for the CASE A and CASE B together. Below 
the curves are the allowed regions. The dashed curves are the 2a allowed region for the 
parameters e and mi. The solid curves are the 3o" allowed region. Thinner curves (black ones) 
represent the CASE A and thicker curves (red ones) represent the CASE B. In Fig. 7 the same 
analysis, now for m4 values. We notice, as shown in both Fig. 6 and Fig. 7, that the value of 
the scale of mi is very similar to m4, since m4 ~ m2, and Am2i ~ ^m\^ ^ 7.5 x 10~^ eV^ 
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for small e. As the radiative correction (e) grows, mi {m^j has to diminish to maintain the 
X^- They have, in some sense, a compensatory behaviour when considered together. We 
remember for the CASE B that e = e,- (e' = e'^), and e and e' are related by Eq. (13). 
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FIG. 6: Allowed regions for the parameters e and mi. Thinner curves (black curves) represent the 
CASE A. Thicker curves (red curves) represent the CASE B. 
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FIG. 7: Allowed regions for the parameters e and m^. Thinner curves (black curves) represent the 
CASE A. Thicker curves (red curves) represent the CASE B. 

In Fig. 6 and Fig. 7, we have plotted all the curves together to make evident the dif- 
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ference between the CASE A and CASE B. We notice that there is a very small difference 
between these two approximations, which would be an evidence that it is almost impossible 
to distinguish between them. 

Our results arc similar to the ones found by de Gouvea et al. [13]. They found in the 
2 + 1 case (two active neutrinos + one sterile) that e < 2.0 x 10"'' for 3a and e < 1.2 x 10~^ 
for 2(7. Their model has only the e parameter, however our model possesses two parameters. 
Despite this fact, as we can see in Fig. 7, we also have found similar values of e when 0.01 eV< 
1714 <0.2 eV at 2a level. The important characteristic of the model, quasi-degenerancy with 
the 1712 state (Am24 much smaller than the other ones), is still mantained even for these 
values of masses. 

V. CONCLUSIONS 

In this paper we have analyzed the model with a quasi-Dirac neutrino put forward in Refs. 
[14, 15] using the solar neutrino data. This is possible because when radiative corrections are 
included in the neutrino mass matrix the oscillation channel Ve — > J^s is open. However, we 
have got that, even in this case, the quasi-Dirac neutrino remains, for all practical proposes, 
a Dirac one, i.e., 1712 ~ m4. Our model has two parameters, the radiative correction e and 
the input mass mi, which is the small one considering, for simplicity, the normal hierarchy. 
We have found allowed regions, shown in Fig. 7, in which e can vary approximately from 
5 X lO-'' to 10-^ and the m4 mass varying from 0.01 eV (0.01 eV) to 0.2 eV (0.3 eV) at 2a 
(3(j) level. These 1714 values are compatible with the most conservative limit of the sum of 
neutrino masses (^^^ < 1.3 eV (95%)) of WMAP-7 [10]. In Ref. [13], which describes a 
2 + 1 model- two active neutrinos plus one sterile - it was obtained e < (1.2,2.0) x 10"^ at 
two and three sigma level, respectively. Note that in our case the four masses belong to the 
interval with order of magnitude ~ (10~^ - 10~^) eV. 

Summarizing, even with radiative corrections are considered, the mixing matrix in the 
lepton sector continues to be the tribimaximal one. It means that in these conditions the 
model cannot explain the disappearance of Pe observed by several experiments and, when 
interpreted in a three active neutrino scenario, it implies a non-zero ^^13 [19-21]. Hence, the 
only way to obtain a realistic PMNS mixing matrix is by considering non-diagonal charged 
lepton mass matrix as has been put forward in Ref. [15]. 
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